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1 .  INTRODUCTION 

In  a  nonequi librium  relativistic  beam-plasma  system  with  a  slowly  varying, 
nearly  spatially  independent  background  distribution  without  external  fields, 
the  total  field,  excited  in  the  beam-plasma  and  associated  with  a  collective 
bremsstrahlung  process  involving  a  test  particle,  has  two  parts.  First,  it 
has  a  regular  part  with  a  relatively  slowly  varying  phase  associated  with  the 
nonradiative  field  of  the  test  particle  and  its  induced  dynamic  polarization, 
and  it  also  has  a  stochastic  part  with  an  irregular  rapidly  varying  phase 
associated  with  the  bremsstrahlung  radiation  field. 1  Thus,  the  Fourier  trans¬ 
form  of  the  total  field  involved  in  the  bremsstrahlung  process  is  given  by 


®k 


♦R 

®k  +  ®k 


(1) 


where  2^  and  2^  are  the  regular  nonradiative  component  and  the  stochastic 
bremsstrahlung  component,  respectively.  The  regular  part,  2£,  of  this  field 
is  assumed  to  consist  of  four  parts  in  sufficient  approximation;  namely, 
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Here  2£* )  is 
relativistic 


the  self -field  associated  with  the  unperturbed  motion  of  the 

- bare  test  particle,  ignoring  its  induced  dynamic  polarization. 

The  field  ^  is  the  field  associated  with  the  perturbation  in  the  bare 
particle  motion.  The  fields  and  are  fields  produced  by  the  dynamic 

polarization  current  induced  by  the  test  particle  to  second  and  third  order  in 
the  total  field,  respectively. 


In  sections  2  through  5,  the  dispersive  characteristics  and  stochastic 
properties  of  the  bremsstrahlung  field  are  reviewed,  the  field  due  to  the 
perturbed  bare  test  particle  trajectory  is  calculated  from  the  associated 
current,  and  the  dynamic  polarization  current  is  used  to  calculate  the  field 
which  it  induces.  In  section  6  an  expression  for  the  total  field  is  obtained, 
and  in  section  7  the  results  are  summarized.  Expressions  for  these  fields  are 
needed  in  calculations  of  the  total  nonlinear  force  on  a  relativistic  test 
particle  involved  in  collective  bremsstrahlung  in  a  nonequilibrium  beam-plasma 


V •  Akopyan  and  V •  N • 
Plasma,  Fiz.  Plazmy,  1  (1975), 


Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
673  [Sov.  J.  Plasma  Phys.,  1_  (1975),  371]. 
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system* The  latter  is  needed  in  calculations  of  collective  bremsstrah- 
lung  processes  and  the  conditions  for  the  occurrence  of  a  bremsstrahlung 
radiative  instability. 


2*  DISPERSION  RELATION  FOR  THE  BREMSSTRAHLUNG  FIELD 


In  zeroth  approximation  the  dispersive  properties  of  the  bremsstrahlung 
radiation  field  sjj*  are  governed  by  the  poles  of  the  linear  photon  Green's 
function  G  (k)  for  the  beam-plasma.  The  latter  determines  the  mode  disper¬ 
sion  relations*  Thus, 


+st  -►o(O)  , 

Efc  *  Efc  /  (3) 

where  the  electromagnetic  waves  of  mode  a  satisfy  the  following  disper¬ 

sion  relation: 


l<Wk>l  -  0  •  (4) 

The  linear  photon  Green's  function  G  (k)  follows  from  the  linear  Maxwell's 

mn 

equations  in  the  beam-plasma  system,  namely, 


V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  1_  (1975),  673  [Sov.  J.  Plasma  Phys.,  1_  (1975),  371]. 

2A.  V.  Akopyan  and  V.  N.  Tsytovich,  Transition  Bremsstrahlung  of  Relati¬ 
vistic  Particles,  Zh .  Eksp.  Teor.  Fiz.,  71_  (1976),  166  [Sov.  Phys.  JETP,  44 
(1976),  87]. 

3A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstralung  Instability  of  Relati¬ 
vistic  Electrons  in  a  Plasma,  Astrofizika,  13  (1977),  717  [Astrophysics ,  13 
(1977),  423]. 

4V.  N.  Tsytovich,  Collective  Effects  in  Bremsstrahlung  of  Fast  Particles 
in  Plasmas,  Comments,  Plasma  Phys.  Conf.  Fusion,  4_  (1978),  73. 

*H.  E.  Brandt,  Nonlinear  Force  on  an  Unpolarized  Relativistic  Test 
Particle  to  Second  Order  in  the  Toted  Field  in  a  Nonequilibrium  Beam-Plasma 
System,  Harry  Diamond  Laboratories ,  HDL-PRL-82-7  (May  1982)  (to  be  published 
as  HDL-TR-1995,  1983). 

*H.  E.  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories , 
HDL-PRL-82-6  (May  1982)  (to  be  published  as  HDL-TR-1994,  1983). 

*H.  E.  Brandt,  Theoretical  Methods  in  the  Calculation  of  the 
Bremsstrahlung  Recoil  Force  in  a  Nonequilibrium  Relativistic  Beam-Plasma 
System,  Harry  Diamond  Laboratories,  HDL-PRL-83-6  (April  1983)  (to  be  published 
as  HDL-TR-2009 ,  1983). 
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The  Fourier  decomposition  of  the  fields  is  given  by 


where 


E  =  /  dk  EkeiO*’?-UtO 
B  -  /dk  , 


dk  =  d#  da)  . 


The  assumed  forms  of  the  constitutive  relations  for  the  beam-plasma  system  are 

+  Bk 

ft,  =  —  ,  (10) 

K  Vo 

where  yQ  is  the  permeability  of  free  space,  and 

°ki  =  eij<*)Ekj  '  <11} 

where  e^j(k)  is  the  dielectric  tensor  for  the  beam-plasma  system,  and  sum¬ 
mation  over  repeated  indices  is  understood.  From  equations  (5),  (10),  and 

(11)  it  follows  that 


\Xq  OcxBk)i  +  iwei j  (^)Ekj  “  3ki 


where  ^  is  the  Fourier  transform  of  the  nonlinear  and  external  current  densi¬ 
ty  ,  and  from  equation  (6) 

B)q  =  — --  ■ 1  —  •  (13) 

u>  +  i6 

Here  6  is  a  small  imaginary  part  of  the  frequency.  Next,  substituting  equa¬ 
tion  (13)  in  equation  (12),  then 

~  [*cx(^xEk)]i  +  iw^ij Ekj  =  jfci  •  (14) 

But  one  has  the  following  vector  identity , 


jfrxftfxEfc)  *  ()t*Ek))?  -  k^Efc  . 


Substituting  equation  (15)  in  equation  (14),  then 
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Polio  +  16) 

Equation  (16)  may  be  revrritten  as 


(k^k^  ~  *2Eki)  +  ^weijEjcj  =  3ki 


(16) 


Gi'j^j 


(0+1 


Z  3ki 


(17) 


where 


,-1 


^  pft«o  +  16)2  ^  i  ) 


(*1*3  -  k2«ij  )  +  tL.  , 


(18) 


where  6.  .  is  the  Kronecker  delta#  For  a  spatially  isotropic  system,  the 
dielectric  permittivity  tensor  is  given  by5'^'7 


6ij(k,io)  =  (fiij  -  ^^et(i?,(0)  +  ejj(k,  (o)  , 


(19) 


where  and  are  the  transverse  and  longitudinal  permittivity,  respec¬ 
tively#  Substituting  equation  (19)  in  equation  (18)  and  combining  terms,  then 


-1  /  k2  \  /  k2  \  ^i^j 

0lj  ■  W  ■  r (TTu)2/Sij  *  VH ' 6* '  fe<«  *  Tst*)  • 


^((U  ♦  i<S)2/  J  >  Wq(< 

Treating  equation  (17)  as  a  matrix  equation,  then 


(20) 


Eki 


to  +  i6 


GijDkj 


(21) 


The  matrix  inverse  to  is  easily  obtained  from  equation  (20);  namely, 

<o2 


ki*j 


'ij 


k*e, 


et(“  +  i«)2  -  IT 


(22) 


Checking  that  equation  (22)  is  in  fact  the  inverse,  one  notes,  using  equations 
(20)  and  (22),  that 


5A.  I.  Akhiezer,  I.  A.  Akhiezer,  R.  V.  Polovin,  A.  G.  Sitenko,  and  K.  N. 
Stepanov ,  Plasma  Electrodynamics ,  _1,  Linear  Theory,  Pergamon  Press  (1975), 
206. 

®V.  N.  Tsytovich,  Nonlinear  Effects  in  Plasma,  Plenum  Press,  New  York 
(1970),  31. 

?V.  N.  Tsytovich,  Theory  of  Turbulent  Plasma  (Consultants  Bureau),  Plenum 
Press,  New  York  (1977),  63-65. 


v  V  H \  r 
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* ■  b-i^Tur^ '  (-  - 


Cki*j 

U*c.  * 


0>* 


l  et(u)  +  i6)2  - 


(•..  -  m  - 


6mj 


as  it  must.  Equation  (21)  can  be  equivalently  written  in  the  following  form: 


***  ’  0)  +  i6  G mn-lkm  ' 


(24) 


since  the  linear  photon  Green's  function  equation  (22),  for  the  isotropic 
plasma  is  a  symmetric  tensor.  According  to  equation  (4),  the  poles  of  this 
function  determine  the  zeroth  order  dispersion  relation  for  the  bremsstrahlung 
in  the  isotropic  case. 

More  generally  the  dielectric  tensor  is  related  to  the  linear  conductivity 
tensor  by 


io^ j  (k) 

e.  .  =  e  6.  .  +  - 

13  0  oj  +  i6 


(25)  i 

i 


Equation  (25)  follows  through  a  comparison  of  equation  (12)  with  the  vacuum 
form  of  Maxwell's  equations;  namely, 


-  (* 

Mo 


)i 


+  ia>eo6ijEkj  =  Oj.  j(k)Ekj  +  jki 


(26) 


where  the  explicit  linear  current  in  terms  of  the  linear  conductivity  has  been 
introduced  on  the  right,*  designates  the  nonlinear  and  external  contribu¬ 
tions  to  the  current,  and  the  vacuum  dielectric  tensor  has  been  introduced. 
Thus,  moving  the  first  term  on  the  right  of  equation  (26)  to  the  left  and 
comparing  with  equation  (12),  then  equation  (25)  follows. 


By  equation  (47)  of  HDL-TR-1994, 
the  beam  plasma  system  is  given  by 


the  linear  conductivity  tensor  a.  .  in 


°ij  =E' 


[*** 
J  (  2ir 


8  vsi 
(28)3 


EU  -  £%)  ♦  M 

- 7 - - T\ - 

i  [(1)  -  k»vs  +  16} 


R(0) 

Ps 


3p 


(27) 


*H.  E.  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories, 
HDL-PRL-82-6  (Nay  1982)  (to  be  published  as  HDL-TR-1994,  1983). 
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Here*  s  is  a  species  label.  Substituting  equation  (27)  in  equation  (25),  then 
the  linear  dielectric  tensor  Cor  the  beam-plasma  system  is  given  by 

,,  ,  _  ,  [*&_  [«,.(’  -  fcgta)  *  ]  “W' 

cij  *  +  A)  +  ifi  *->e®  )  (2v)  3  (w  -  k»vs  +ifi)  ^Psm  '  (28) 


in  which  case  equation  (18)  becomes 


37]  * - 5 - rr  fak*  -  k2^  )  +  tn 6,, 

u  (to)  +  i«)2  v  1  ^  i3;  °  x3 


+  •"i"i*?es, 


vsi 


r,  (, .  *3l_)  *  M*<o) 

L3m^  o)  +  at  +  jgJ  P3 


(1)  -  +  i6 


It  should  be  noted  that  equations  (27),  (28),  and  (29)  are  in  fact  symmetric 
in  i  and  j  which  follows  after  a  simple  integration  by  parts.  The  zeros  of 
the  determinant  of  the  matrix  equation  (29)  must  give  the  dispersive  prop¬ 
erties  of  the  bremsstrahlung  field. 

From  equation  (19)  one  has  that 


k2  6ij  ' 


Et  *  1  (eii  "  ef)  * 


Substituting  equation  (28)  in  equations  (30)  and  (31),  then  in  the  spatially 
isotropic  case,  the  longitudinal  and  transverse  dielectric  permittivity  are 
given  by 


1  f  d3PS  k*Vs  ,♦  ♦ 

®i  *  eo  +  - I - 0  - * -  k«V, 

u>  +  16  s  k2  J  (2ir)3  at  -  k»vs  +  ifi 


;  fR(o) 

PsfPs 


®t  -  ®o 


1  f  d3Ps  T  o  k2v|  -  oic  «v3  ~|  3fPg( 

2 (to  +  ifi)2  a  k2  J  (2»)3  L  8m  10  -  )c*Vfl  +  ifi  ^J^Psm 


If  the  background  distribution  functions  f£(°)  are  also  isotropic  in  momentum 
space,  then  they  depend  only  on  particle  energy;  namely, 

fS(0>  "  fS(0>(O  »  (34) 


10 


where 
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(35) 


Then 
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Substituting  equation  (36)  in  equation  (28)  and  simplifying ,  then  for  the  case 
of  isotropy  in  momentum  space,  one  has 


Hi  -  — ^£-1  [&.*?)  - 

o)+i6~  y  ( 2w )  (u  -  k»vg  +  1 


3f 


R(0) 


Ps 


iS) 


3ec 


(37) 


Substituting  equation  (37)  in  equations  (30)  and  (31),  then  in  the  case  of 
isotropy  in  momentum  space,  the  longitudinal  and  transverse  dielectric  permit¬ 
tivity  are  given  by 


e* 


. .  ,  i  Yd  fil* 

®  id  +  i6^-*  k2  J  (2n 


S-v.)s 


3f 


R(0) 


Ps 


(2w)3  u  -  k *vs  +  i6 


(38) 


and 


1  '  v-  ,  fd3Ps  *1  -  (it*va)2/k2  <0> 

Ct“°*tzrnZeiJ  Tsn 


(2tt)3  0)  -  k«vs  +  i6  3es 


(39) 


and  the  zeros  of  the  determinant  of  equation  (20)  along  with  equations  (38) 
and  (39)  then  give  the  dispersion  relation  for  the  bremsstrahlung  in  this 
case* 


THE  STOCHASTIC  PROPERTIES  OF  THE  BREMSSTRAHLUNG  FIELD 


As  stated  earlier,  it  is  to  be  assumed  that  the  bremsstrahlung  field  re¬ 
sulting  from  the  particle  scattering  in  the  system  is  stochastic,  having  an 
irregular  rapidly  varying  phase.1  To  the  needed  order,  the  phase  of  the 
stochastic  bremsstrahlung  field  is  assumed  to  change  randomly  during  a  typical 
interaction  period,  and  the  time  average  is  assumed  to  be  the  same  as  the 
average  over  the  stochastic  ensemble.  In  particular  it  is  assumed  that  for 
this  stochastic  field  the  following  averages  over  the  statistical  phase 
distribution  apply  approximately;  namely. 


1A.  V •  Akopyan  and  V •  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazmy,  1_  (1975),  673  [ Sov .  J.  Plasma  Phys.,  J_  (1975),  371]. 
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where  e^  is  the  unit  polarization  vector  for  the  stochastic  field. 


Equations  (40)  and  (41 )  are  deduced  as  follows.  Representing  the  stochas¬ 
tic  electric  field  in  terms  of  its  modulus,  A^,  and  phase,  one  has 


,st 


®ki  “  ekiAke 


... 


as  e 


kiAk(cos 


+  i  sin 


*k) 


(42) 


Therefore,  the  average  over  the  statistical  phase  distribution  <E^^>  is  given 
by 


<Eki>  “  ekiAk (<cos  *k  +  i  sin  *k>)  • 

But  assuming  random  phase  to  the  needed  order,  then  clearly 


(43) 


<cos  4)c>  =  <sin  (^>  =  0 


(44) 


and  equation  (43)  becomes  equation  (40).  Similarly, 


at-  a*.  ”i  ft,  i 

<EkiEk'j>  =  eki*jc,jAkAk,<e  e  > 


(45) 


But 


<.lV*V 


>  ^ 


<e  *  K  >  •  6^, 


(46) 


since  it  is  vanishing  unless  k  =  k*  and  unity  when  k  =  k'.  Substituting 
equation  (46)  in  equation  (45),  using  the  property  of  the  Kronec^cer  delta, 
denoting  k*  by  and  using  the  reality  property  of  the  field  =  E^), 

then 


<EkiEk^j>  ^i^j^k^-k^  ,k  * 

Next,  one  may  define  the  quantity  |E^t|  in  terms  of  the  modulus  A^  by 

l*k*|2«(*  .  (48) 

Then  substituting  equation  (48)  in  equation  (47),  equation  (41)  follows. 


Furthermore,  using  equation  (3)  in  equations  (40)  and  (41)  one  has  for  the 
stochastic  properties  of  the  bremsstrahlung  radiation  field  the  following: 


<Eki0)>  =  0 


(49) 


<Ek°i°)<?)>  =  |Ek°(0)|2  6(k  +  k,  ) 


where  e£  Is  the  unit  polarization  vector  for  mode  o. 


(50) 


4*  THE  FIELD  DUE  TO  THE  PERTURBED  TEST  PARTICLE  TRAJECTORY 

The  trajectory  of  a  relativistic  test  particle  of  species  a  moving  through 
the  beam-plasma  system  is  given  by 

ra  =  vat  +  Ara  ,  (51 ) 

where  vfl  is  the  unperturbed  velocity  and  and  Av^  describe  the  perturbation 
in  the  trajectory  and  the  velocity,  respectively.  Provided  the  conditions  for 
the  plasma  Born  approximation  are  satisfied — namely,  that  the  relativistic 
particle  momentum  be  much  greater  than  the  electromagnetic  impulse  received  by 
the  particle  in  a  time  interval  given  by  the  inverse  electron  plasma  frequen¬ 
cy — then  the  perturbed  motion  is  given  by  equation  (36)  of  HDL-TR-1 995: * 


n  , 
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-i[u>-k*v  Jt 

va#Ek  1  e 
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w  +  i6Ja)-k»v  +16 


The  current  density  due  to  the  actual  charge  of  the  test  particle  is  then 
given  by 


=  ea(va  +  Ava)63(r  -  vat  -  Ara) 


Expanding  the  delta  function  in  equation  (56)  for  small  A?a,  then  to  lowest 

_ j _ j  _ 


order  in  ArQ, 


la  =  ea^a  +  A*a)[63(*  ”  V)  '  ^a*763^  '  V)1 


Thus,  to  lowest  order  in  A?a  and  Ava,  one  has 


*  .  *<1>  +  +(2) 

Ja  Ja  Ja 


where 


n"  ■  -  V) 


Da2)  =  eaAva63(r  -  vat)  -  eavaAra»$53(r  -  vat) 


The  Fourier  components  of  these  currents  at  the  position  of  the  test  particle 
are  given  by 


i&l  - / 

)£l-f • 


Substituting  equation  (59)  in  equation  (61  )  and  integrating  over  the  delta 
function,  then 


gYcd  f  dt  -i(^*v0(-<0)t 
2n)3j1" 


Next,  recognizing  the  integral  form  of  the  delta  function  appearing  in  equa¬ 
tion  ( 63 ) ,  then 
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Next,  substituting  equation  (60)  in  equation  (62)  one  has 


-  /  eavod^a*[^3(^  -  vat)  ]e-i  (k*r-«t)  . 


(65) 


The  second  integral  in  equation  (65)  can  be  simplified  by  integration  over 
space  by  parts  and  dropping  surface  terms.  Equation  (65)  becomes 
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Integrating  over  the  delta  functions,  then 
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Next,  substituting  equations  (52)  and  (55)  in  equation  (67)  and  integrating 
over  the  time,  then 


.(2)  /vcivaj  „  \  /  d^du,  [_  (  *i *va  ^ 

Ddki  c2  “  6ij/  i  (2tt)  3  LEkin  "  “1  ♦  W 

~|  «(«  -  m,  -  (*  -  1?,  HJ 

+  1j  u)1  +  U]  .  f  +  +  i6 

1  i  a  jt 


ie2  f  d3)^  day,  £**^k,0  ay,  +  ii)  + 

+  V(xiJ  ~72^  "(«,  -  Jt1.^a+  ay 


(68) 


0). 


”T5'J 


x  6(a>  -  (it  -  ic1  )*vj 


ie2 

a 


y  cPcF 


2  ai 


f  dtf.du).  *•*«  v^. 

I  - — - - — r - - - —  6  fw  -  w-  -  fit  -  it.  )*vrt)  • 

J  (2ir)3  (a),  -  <t,4a+  ifi)2  1  11  lJ  aJ 


Simplifying  equation  (68),  then 
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oj1  -  k1  •va  +  16  -  k1  «va  +  i6 Jz  c^[w1  -  k^  •v(J  +  i 6 J ^ 

Usinc^  the  property  of  the  delta  function  to  replace  (u)  -  k  «v  ) 
(w  -  £*va)  in  the  fourth,  fifth,  and  sixth  terms  of  equation  (69),  then  a 
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Changing  the  integration  variable  k1  =  ( 0>1 , )c  1  )  to  -k1 ,  equation  (70)  becomes 
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Combining  terms  and  using  the  property  of  the  delta  function,  one  notes  that 


VV^FJ '  •  v°~  ’vv°t  • 

Substituting  equation  (72)  in  equation  (71)  and  using  the  reality  property  of 
the  field  (e^  *  then 


(73) 
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where 

The  currents  and  given  by  equations  (64)  and  (73)  are  the  cur¬ 

rents  associated  with  the  unperturbed  motion  and  the  perturbation  in  the 
motion,  respectively,  of  the  bare  particle.  By  the  same  arguments  leading  to 
equation  (21),  the  components  of  the  field  due  to  the  unperturbed  and  pertur¬ 
bed  test  particle  trajectory  are  given  by 


)6ij  + 


vcdkj  “  k1ivaj 


UXl)- 

vniv  •  (*•*! - r) 

OO-  Of)  »■  1  C2J 


(74) 
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j?«v0  -  i6  (a)  -  it*va  -  i6)2 


p(  1 )  _ 
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oj  +  i6 


(1) 


'run 


3km 


(75) 


and 


e>( 2)  _ 
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u>  +  i( 


Gnm3km^ 


(76) 


respectively.  Therefore,  substituting  equations  (64)  and  (73)  in  equations 
(75)  and  (76),  one  obtains 


ie. 


(2n)3(u>  +  i6)  m  nm 


(k ) 6(u 


(77) 


and 


*ln 


- -  G  (k ) 

(2tt)3(o)  +  i  6)  nm 


)t*v  -  •*  ) 

a  i  • 
(78) 


Equations  (77)  and  (78)  agree  with  equations  (16)  and  (17)  of  Akopyan  and 
Tsytovich,  1975,  1  except  for  erroneous  overall  factors  of  ( 2ir) ~ 3  and  (2n)~3i 
appearing  there  in  equations  (16)  and  (17),  respectively.  The  indices  n  and  m 
are  interchanged  there  as  they  are  in  the  defining  relation,  equation  (14)  of 
Akopyan,1  compared  to  equation  (21)  here.  However,  a  spatially  isotropic 
plasma  is  assumed  there,  so  according  to  equation  (22)  the  Green's  function  is 
symmetric,  namely  G,^  =  Gnm.  Also  in  Akopyan,  1  single-wave  particle  resonance 
is  ignored,  in  which  case  the  small  imaginary  part  iS  in  equation  (74)  is 


*A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsst rah  lung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plazm y,  1  (1975),  673  [Sov.  J.  Plasma  Phys.,  1  (1975),  371], 
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/ 

t,  iLj  becomes  effectively  real,  and  the  complex  conjugate  sign  in 
(78)  may  be  removed.  The  additional  factors  of  (4s) (2s) 3  in  equation 


ignorable, 
equation 

(16)  and  (4s)  in  (17)  there  are  due  to  the  Gaussian  form  of  Maxwell's  equation 
used  there  compared  to  the  MKS  form  used  here,  and  to  the  difference  in  Fouri¬ 
er  transform  convention.  For  example,  in  equation  (5)  of  Akopyan  and 
Tsytovich,  1975, *  the  Fourier  transform  convention  employed  has  a  factor  of 
(2s)  ^  in  the  inverse  Fourier  transform  in  the  integration  over  the  three- 
dimensional  wave  vector  space,  and  a  factor  of  1  for  the  integration  over 
frequency,  giving  a  total  factor  of  (2s)'3.  In  this  report,  a  total  factor  of 
1  is  used,  as  in  equation  (7),  for  example.  Likewise,  the  Fourier  transform 
itself  then  has  a  factor  of  (2s)'1  there  and  (2s)'4  here,  as  in  equation  (61), 
for  example. 


5.  THE  FIELD  DUE  TO  THE  DYNAMIC  POLARIZATION  CURRENT 

The  Fourier  transform  of  the  dynamic  polarization  current  to  third-order 
induced  in  the  beam-plasma  is  given  by1' 


t  =  ?(1  >  +  t<2)  +  *(3) 

3dpk  3dpk  3dpk  +  3dpk 

The  linear  polarization  current  density  J^pk  ven  by 


3dpki  “  I  °i!)  WEk* 


(s) 


(79) 


(80) 


Its  effects  are  manifest  in  the  linear  Green's  function  for  the  beam-plasma 
system  as  discussed  in  section  2.  The  second-order  and  third-order  polar¬ 
ization  current  densities  and  are  given  by 


ji22  =  -T,e  ffH  dk26^  ~  *1  ~  hi  =(s )fl  k  .  i 

3dpki  4-  s  J  +  +  i6)  ' 


(81  ) 


and 
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(82) 


xZijfm(k'k1'k2'k3)Ek1jEk2fEk3m  ' 

respectively*  The  second-order  nonlinear  and  third-order  nonlinear  conductiv- 
ities  for  species  s,  and  E^®^,  respectively,  are  given  by 


1A.  V.  Akopyan  and  V .  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fix.  Plaza y,  (1975),  673  [ Sov .  J.  Plasma  Phys.,  1_  (1975),  371], 

*H.  E*  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories, 
HDL-PRL-82-6  (May  1982)  (to  be  published  as  HDL-TR-1994,  1983). 
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The  function  fj^°)  is  the  assumed  slowly  varying  and  spatially  independent 
background  distribution  function  for  species  s,  whose  Fourier  transform  is 

fk(S)  =  f£(0)5<k>  •  (85) 

K  ps 

To  determine  the  fields  defined  in  equation  (2)  and  associated  with  the 
respective  parts  of  the  dynamic  polarization  current,  one  proceeds  as  in 
obtaining  equations  (75)  and  (76).  Thus,  one  has  for  the  fields  produced  by 
the  second-  and  third-order  nonlinear  dynamic  polarization  currents,  respec¬ 
tively. 
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Next,  substituting  equation  (81)  in  equation  (86),  one  obtains 
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First  symmetrizing  equation  (88)  in  ( j )  and  (t,k2),  changing  variables  of 
integration  from  k^  to  -k1  and  k2  to  -k2,  and  using  the  reality  property  of 
the  field,  then  equation  (88)  becomes 


B^n  “  2( w  +  i«>  Gn»(k)  ?*8  /  («,'  -  i«)(«22-  TT)  6(k  +  kl  +  k^ 
x  '"k2)  +  Smlj  )JEk1jEk2£  • 


(89) 


Comparing  equation  (89)  with  equation  (22)  of  Akopyan,1  a  disparity  of  a 
factor  of  (1/2)  (4ir)  (2ir)“^  in  the  latter  j.s  apparently  due  to  explicit  inclu¬ 
sion  there  of  a  self -field  factor  of  (1/2),  the  use  of  Gaussian  units,  and  the 
different  Fourier  transform  and  normalization  conventions  chosen  there.  The 
factor  of  (1/2)  is  either  a  typographical  error  or  results  from  explicit 
inclusion  of  a  factor  of  1/2  associated  with  self-fields.  For  the  generic 
expression,  equation  (89),  it  is  preferable  that  such  a  factor  be  maintained 
implicit  until  the  formula  is  applied  explicitly  and  a  self -field  factor  of 
(1/2)  is  needed.  The  factor  of  (4tt)  is  due  to  the  use  of  Gaussian  units 
there.  One  of  the  factors  of  (2ir)“3  is  due  to  the  differing  Fourier  transform 
convention.  Because  of  the  different  Fourier  transform  convention,  the  coun¬ 
terpart  of  equation  (81)  would  have  an  additional  factor  of  (2tt)“3.  The  other 
factor  of  (2ir)”3  is  apparently  due  to  the  different  normalization 
of  fpg0^  there.  In  short,  the  fp^0^  there  is  evidently  ( 2 ir )  ^  times  that 
here.  Alternatively,  if  the  normalization  is  the  same  as  that  here  then  there 
is  an  erroneous  factor  of  (2x)~3  appearing  there.  An  isotropic  plasma  is 
assumed  there,  in  Which  case  by  equation  (22)  the  Green's  function  is  symmet¬ 
ric,  namely,  =  Gnm.  However  as  already  noted  the  indices  of  G^ .  are 

interchanged  in  Akopyan  and  Tsytovich1  in  the  definition  of  the  Green's  func¬ 
tion.  It  appears,  however,  that  there  are  typographical  errors  there.  The 
complex  conjugate  sign  on  the  fields  and  a  factor  of  eg  are  left  out.  Also, 
there  in  the  second-order  nonlinear  conductivity,  one  notes  from  equation  (83) 
that* 


smj  tO^-kl  '**2)  =  smjt  (-*#*1**2)  =  ~^mjf  (k,kj,k2) 


,(s) 

*mj 


«<s)< 

mj  J 
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where  ^(k,k-|  ,kj)  designates  the  conductivity  tensor  appearing  in  Akopyan 

and  Tsytovich.  In  the  last  step  of  equation  (90)  it  is  recognized  that 

3^^(k,ki,k2)  differs  from  smj  ,k2)  here  in  that  the  first"  complex  de¬ 
nominator  <0  -  £*v  +  i6  in  equation  (83)  here  is  implicitly  to  -  -  i6 


*A.  V .  Akopyan  and  V .  N .  Tsytovich ,  Bremsstrahlung  in  a  Nonequilibrium 
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there. 8# 9**  Thus,  there  is  also  a  disagreement  in  overall  sign  between  equa¬ 
tion  (89)  here  and  equation  (22)  of  Akopyan  and  Tsytovich.  1  In  summary 
then,  equation  (22)  of  Akopyan  and  Tsytovich1  is  in  error;  it  should  have  an 
additional  factor  of  (-eg),  the  fields  and  the  conductivity  tensor  on  the 
right-hand  side  should  be  complex  conjugated,  and  the  explicit  additional 
factor  of  (1/2)  there  should  be  omitted. 


Next,  substituting  equation  (82)  in  equation  (87), 
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Symmetrizing  equation  (91)  in  (i,k*,),  (j,k2)>  (f,k3),  changing  variables  of 

integration  from  k^  to  -k^,  k~  to  -k2,  and  k^  to  -k^,  and  using  the  reality 
property  of  the  field  ^  J,  then  equation  (91)  becomes 


(92) 


Comparing  equation  (92)  with  equation  (23)  of  Akopyan  and  Tsytovich,  1  the 
disparity  of  a  factor  of  ( 1/2)  (4tt)  (2tt)"9  in  the  latter  is  again  due  to  an 
explicit  self-field  factor  of  (1/2),  the  use  of  Gaussian  units,  and  the  dif¬ 
ferent  Fourier  transform  and  background  normalization  conventions  there. 
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Alternatively  if  the  normalization  is  instead  the  same  as  that  here,  then 
there  is  an  erroneous  factor  of  (2tt)“3  appearing  there.  There  is  apparently 
an  error  there.  A  factor  of  ies  has  been  left  out.  Note  that  by  equation 
(84)  it  follows  that 


Z  mijA^k,_k1  '”k2'-k3  ^  EmijJl(-k'k1  ,k2'k3) 
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(93) 


conductivity  tensor  T$?  '^2*^3 )  9iven  bY  equation  (21)  of 

Tsytovich,  1  the  first  two  complex  denominators  are  implicitly  a 
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Also  in  the 
Akopyan  and  Tsytovich 
-  i?*vs  -  i6  and  oj  +  031  -  ()c  +  )  *v^  -  i6,  respectively,  there.*  Also,  there 

is  no  overall  factor  of  i  as  there  is  in  equation  (84)  here.  Therefore, 
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Combining  equations  (93)  and  (94)  then 
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In  summary  then,  equation  (23)  of  Akopyan  and  Tsytovich1  is  also  in  error  and 
should  also  have  an  additional  factor  of  (ieg),  the  complex  conjugate  of  the 
fourth-order  conductivity  tensor  there  must  be  taken,  and  the  explicit  addi¬ 
tional  factor  of  (1/2)  there  should  be  omitted. 


6.  THE  TOTAL  FIELD 


Collecting  equations  (1),  (2),  and  (3)  the  total  field  involved  in  the 
bremsstrahlung  scattering  process  and  acting  on  the  test  particle  is  given  by 


K  -  K'°'  +  *  K.2)  +  ■ 


(96) 


Here 


o(0) 


is  the  bremsstrahlung  field  whose  dispersion  relations  are  given  by 


equations  (4)  and  (29)  in  the  spatially  isotropic  case  and  equations  (4), 
(20),  (38),  and  (39)  in  the  case  of  isotropy  in  both  ordinary  space  and  mo¬ 

mentum  space.  The  stochastic  properties  are  given  by  equations  (49)  and 
(50).  The  fields  iS^1 '  and  given  by  equations  (77)  and  (78)  are  the 

self -fields  arising  from  the  relativistic  test  particle's  own  motion.  The 
fields  *dpk and  *dpk  given  by  equations  (89)  and  (91)  are  the  increasing  order 
fields  produced  by  the  nonlinear  dynamic  polarization  current  induced  by  the 
test  particle.  The  fields  and  are  excited  in  nonlinear  scattering 

of  the  field  by  the  induced  aynamic  polarization  current.  They  involve  the 
three-  and  four-plasmon  vertex,  respectively,  and  play  a  key  role  in  deter¬ 
mining  the  nonlinear  contribution  to  the  collective  bremsstrahlung  process. 


*A.  V.  Akopyan  and  V .  N.  Tsytovich ,  Bremsstrahlung  in  a  Nonequilibribrn 
Plasma ,  Fiz.  Plazmy ,  1_  (1975),  673  [ Sov .  G.  Plasma  Phys.,  _2_  (1975),  371]. 

*H.  E.  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories, 
HDL-PRL-82-6  (May  1982)  (to  be  published  as  HDL-TR-1994 ,  1983). 
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7*  CONCLUSION 


An  expression  has  been  obtained,  equations  (96),  (77),  (78),  (89),  and 
(91),  for  the  total  field  involved  in  the  collective  scattering  and  brems- 
strahlung  of  a  relativistic  test  particle  in  a  nonequilibrium  beam-plasma 
system*  Dispersion  relations,  equations  (4),  (29),  (20),  (38),  and  (39),  and 
stochastic  properties,  equations  (49)  and  (50),  for  the  bremsstrahlung  field 
have  also  been  obtained*  These  results  are  in  agreement  with  those  of  Akopyan 
and  Tsytovich1  with  the  exception  of  apparent  typographical  errors,  and  have 
been  used  by  them  in  their  work  on  nonlinear  bremsstrahlung  in  nonequilibrium 
plasmas* 

The  present  work,  together  with  related  work  by  the  author,  *++§  is 
important  for  ongoing  work  in  calculating  collective  radiation  processes  and 
conditions  for  the  occurrence  of  radiative  instability  in  relativistic  non¬ 
equilibrium  beam-plasma  systems. 


lA.  V.  Akopyan  and  V •  N.  Tsytovich ,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz .  Plazmy,  (1975),  673  [ Sov .  J.  Plasma  Phys.,  1_  (1975),  371]. 

10tf*  E.  Brandt,  Symmetries  of  the  Nonlinear  Conductivity  for  a  Relativistic 
Turbulent  Plasma,  Harry  Diamond  Laboratories,  HDL-TR-1927  (March  1981). 

1  ^H.  E.  Brandt,  Exact  Symmetry  of  the  Second-Order  Nonlinear  Conductivity 
for  a  Relativistic  Turbulent  Plasma,  Phys.  Fluids,  24_  (1981),  1760. 

12H*  E.  Brandt,  Second-Order  Nonlinear  Conductivity  Tensor  for  an  Unmagne¬ 
tized  Relativistic  Turbulent  Plasma,  in  Plasma  Astrophysics,  Course  and 
Workshop  Organized  l*y  the  International  School  of  Plasma  Physics,  27  August  to 
7  September  1981,  Varenna  (Como),  Italy  (European  Space  Agency  ESA  SP-161, 
November  1981),  361  (also  to  be  published  by  Pergamon  Press). 

13tf*  E.  Brandt,  Symmetry  of  the  Complete  Second-Order  Nonlinear  Conductiv¬ 
ity  Tensor  for  an  Unmagnetized  Relativistic  Turbulent  Plasma,  Journal  of  Math. 
Phys .  2±,  (1983),  1332,  2250. 

^H.  E.  Brandt,  The  Gluckstern-Hull  Formula  for  Electron-Nucleus  Brems¬ 
strahlung,  Harry  Diamond  Laboratories,  HDL-TR-1884  (May  1980). 

*H.  E.  Brandt,  Nonlinear  Force  on  an  Unpolarized  Relativistic  Test  Parti¬ 
cle  to  Second  Order  in  the  Total  Field  in  a  Nonequilibrium  Beam-Plasma  System, 
Harry  Diamond  Laboratories ,  HDL-PRL-82-7  (May  1982)  (to  be  published  as  HDL- 
TR-1995,  1983). 

t H .  E.  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories, 
HDL-PRL-82-6  (May  1982)  (to  be  published  as  HDL-TR-1994,  1983). 

Tff*  E.  Brandt,  Theoretical  Methods  in  the  Calculation  of  the 
Bremsstrahlung  Recoil  Force  in  a  Nonequilibrium  Relativistic  Beam-Plasma 
System,  Harry  Diamond  Laboratories,  HDL-PRL-83-6  (April  1983)  (to  be  published 
as  HDL-TR-2009,  1983). 

* Other  related  work  prepared  in  preprint  form  will  be  published  later  and 
is  available  from  the  author. 
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